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Are cosmological neutrinos free-streaming? 
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Precision data from cosmology suggest neutrinos stream freely and hence interact very weakly 
around the epoch of recombination. We study this issue in a simple framework where neutrinos 
recouple instantaneously and stop streaming freely at a redshift Zj. The latest cosmological data 
imply Zi < f500, the exact constraint depending somewhat on the assumed prior on z%. This bound 
translates into a limit on the coupling strength between neutrinos and majoron-like particles <f>, 
implying r > f x f 10 s (m2/50 meV) 3 for the decay V2 —* v\ + <f>. 

PACS numbers: 98.80.-k, 14.60.St, 14.80.Mz 



I. INTRODUCTION 

With the advent of high-precision cosmology it has be- 
come feasible to probe progressively more detailed as- 
pects of the cosmic neutrino background radiation [l], Q ■ 
In the standard model, neutrinos provide relativistic en- 
ergy density which influences the cosmic microwave back- 
ground (CMB) radiation mainly via the early Integrated 
Sachs-Wolfe (ISW) effect and the matter fluctuation 
spectrum via the relation between neutrino energy den- 
sity and the epoch of matter-radiation equality. The 
existence of a cosmological background of relativistic en- 
ergy density has been unambiguously detected in the 
WMAP-5 data [3J and was already previously detected 
using the combination of CMB and Large Scale Structure 
(LSS) data [I HIHl. Furthermore, cosmological data 
provide a restrictive upper bound on the sum of neutrino 
masses of 0.2-1 eV, depending on the specific choice of 
data sets and model space 0, [J G3, M, Qj, [H, El El El] • 

The present level of precision allows us to turn to 
more subtle issues. For example, it is timely to probe 
the possibility that neutrinos have non-standard interac- 
tions where one case in point is an interaction with the 
Nambu-Goldstone boson of a new, broken 17(1) symme- 
try as in majoron models [l?], El Ell- Such an interac- 
tion would recouple the neutrinos to each other at some 
"interaction redshift" Zi, whereas at earlier epochs they 
would behave in the same way as standard-model neu- 
trinos. For the cases of interest, this recoupling occurs 
much later than the electroweak decoupling. Therefore, 
in the limit of relativistic neutrinos the total energy den- 
sity in the combined fluid of neutrinos and majorons is 
conserved, preventing any direct impact on cosmological 
observables. 

However, neutrinos lose their free-streaming property 
if the interaction is sufficiently strong. As a consequence, 
any anisotropic stress components in the Boltzmann hi- 
erarchy are suppressed, effectively truncating the Boltz- 
mann hierarchy at first order, equivalent to the equations 
for a perfect fluid [10, HH HI, H3, H3, HE HI, H3, Hi] . (See 
Ref. [2!| for a detailed description of the Boltzmann hi- 
erarchy.) 



The impact of neutrino free streaming on cosmolog- 
ical observables was recently studied in Ref. 28]. The 
fit parameter was the effective viscosity c^ is , taken to 
be independent of redshift. Our study is complemen- 
tary in that we assume that c^ is drops instantaneously at 
Zi from the free-streaming value 1/3 to the perfect-fluid 
value OL Our conclusion that neutrinos should stream 
freely around the epoch of recombination is perfectly con- 
sistent with Ref. [2J|. However, our approach lends itself 
more directly to an interpretation in terms of a specific 
interaction model where the recoupling redshift is related 
to a dimcnsionless coupling constant g. Therefore, we can 
translate our limits on Zi into limits on g. 

From flavor oscillation experiments we know that neu- 
trinos have masses which therefore are unavoidable cos- 
mological fit parameters. The usual cosmological lim- 
its on the sum of neutrino masses imply that any single 
mass eigenstate should obey m < 0.2-0.3 eV so that 
all neutrinos would be relativistic around the recombina- 
tion epoch. Treating them as massless is therefore a rea- 
sonable approximation for the simple problem addressed 
here. On the other hand, a strong majoron-type inter- 
action can lead to the annihilation of "heavy" neutrinos 
into majorons ( "neutrinoless universe" (22J). Such sce- 
narios lead to a complicated evolution of the neutrino- 
majoron fluid that we are not investigating, although it 
would have a strong impact on cosmological observables. 
In any event, our constraint on the free-streaming na- 
ture of the relevant radiation at recombination does not 
depend on the physical nature of the radiation. 

Eventually the KATRIN experiment, unless it detects 
a significant neutrino mass, will constrain the neutrino 
mass scale to m < 0.2 eV 30]. Such a bound would 
imply that neutrinos can not have disappeared at the 
recombination epoch and our constraint indeed applies to 
neutrinos. In this sense the anticipated KATRIN limit 
will strengthen the case for translating our limit on Zi 
into a limit on exotic neutrino interactions. 



t A paper by the same authors with a more extensive treatment 
of this issue is in preparation. 



2 



We begin in Sec. [TT] with a description of our model 
space, data sets, and statistical methodology. In Sec. ITTT1 
we provide our bounds on Zi that are translated, in 
Sec. IIV1 into limits on a neutrino-majoron coupling 
strength g. We conclude in Sec. IVl 



II. MODELS, DATA, AND METHODOLOGY 

Our parameter constraints will be based on a reason- 
ably general 10-parameter model consisting of 

9 = (u>ct>m,ub, H ,n s ,a s ,T, A s , N v , z^, (1) 

where h = iJo/(100 km s^ 1 Mpc -1 ), and the cold 
dark matter and baryon contents are given by wcdm = 
JIcdm^ 2 and wb = fie ft 2 respectively. We assume spa- 
tial flatness, i.e. the dark energy density is given by 
^de = 1 — Ocdm — ■ For the dark energy we assume a 
constant equation of state parameter w. The primordial 
fluctuations are assumed to be adiabatic and described 
by the scalar amplitude A s , the spectral index n s , and 
the running a s . We do not consider the presence of tensor 
modes or an isocurvature component. Finally, we include 
the present Hubble parameter, Hq, and the optical depth 
to reionization, r. As discussed in the previous section 
we assume massless neutrinos. 

In order to keep our study on the properties of the ra- 
diation as general as possible, we will sometimes use the 
effective number of neutrino flavors, N v , as a fit param- 
eter to express the radiation content in the usual way. 
The standard value is N v = 3.046 [S3. 

Neutrino interactions are assumed to recouple instan- 
taneously at a redshift Zj. Here our standard prior is 
linear (i.e. uniform) in z^, but we will also test an alter- 
native case where a linear prior is used on log(^ + 1). 

The priors on our model parameters are listed in Ta- 
ble|Tl including the alternatives that we use in some cases. 



TABLE I: Priors on the model parameters. Note that for 
N v we use the prior < N v < 50 for Figs. [1] and [2] and 
N v = 3.046 everywhere else. 

Parameter Standard (linear) prior Alternative (log) prior 



WCDM 


0.01-0.9 


0.01-0.9 




0.005-0.1 


0.005-0.1 


Ho 


40-100 


40-100 


w 


-2-0 


-2-0 


n s 


0.5-1.5 


0.5-1.5 


a a 


-0.2-0.2 


-0.2-0.2 


logllO 10 ^] 


2.5-4 


2.5-4 


T 


0-1 


0-1 


N v 


3.046 / 0-50 


3.046 / 0-50 


Zi 


0-10 4 




log(l + Zi) 




0-4 



We use CMB data from WMAP-5 0, IH H, S| and 
measurements of the matter power spectrum based on 
the Sloan Digital Sky Survey-Luminous Redshift Galax- 
ies (SDSS-LRG) gf| and 2-degree-Field (2dF) galaxy 
samples [35| . In addition we include the Supernova 
Type la (SN-Ia) data from Ref. [H, the SDSS-LRG 
Baryon Acoustic Oscillation measurement (SDSS-LRG 
BAO) from Ref. [H], and the Hubble Space Telescope 
(HST) key project measurement of H pfJ |. 

Our treatment of non-linear corrections to the LSS 
power spectra follows the prescription given in Ref. [4l| . 
i.e., we include data up to k = 0.2/iMpc -1 and correct 
for non-linearity using the shot-noise term P s hot- 

In order to derive constraints on our model param- 
eters we have modified the publicly available CAMB 
code p5l [ to allow for neutrino interactions and com- 
bined it with the Markov Chain Monte Carlo software 
COSMOMC Credible intervals are calculated using 
Bayesian inference as implemented in the GetDist routine 
of COSMOMC. 



III. LIMIT ON RECOUPLING REDSHIFT 

Following the approach described in the previous sec- 
tion we have calculated 68% and 95% credible regions 
in the 2D parameter space of N v and Zi that we show in 
Fig.[TJ In the upper panel we have used the linear prior on 
Zi described in Table Q] The cosmological precision data 
show (i) strong evidence for the existence of relativistic 
energy density and (ii) that it must be freely streaming at 
a redshift around recombination (z r 1100). Marginal- 
izing over N v we find Zi < 1500 at 95% C.L. 

We have repeated the same exercise for a logarithmic 
prior on z,, i.e., one that is uniform in log(zi + 1). The 
corresponding credible regions are shown in the lower 
panel of Fig. [1] Marginalizing once more over N v we 
find Zi < 795 at 95% C.L. The difference arises because 
the effective volume at low Zi becomes larger for the log- 
arithmic prior and therefore integration favors slightly 
lower values of Zi . This effect is well known for parame- 
ters with a highly non-Gaussian likelihood, other notable 
examples being neutrino mass, rriy [lq ], and the tensor 
to scalar ratio, r [43ll44l|. 

While the exact redshift at which neutrinos can be- 
come strongly interacting depends on assumptions about 
priors, we find that neutrinos which were strongly inter- 
acting significantly before recombination are excluded by 
data at much more than 95% C.L., a conclusion which is 
fully consistent with Ref. f28j . 

Our results pertain to any form of radiation present 
around recombination. However, we ultimately want to 
test the interactions of ordinary neutrinos. The cosmic 
standard radiation content is given by N v — 3.046. Re- 
peating the above exercises with this fixed prior we find 
Zi < 1520 for the linear Zi prior and Zi < 790 for the log- 
arithmic prior. These limits are almost identical to those 
where we marginalized over N v . This is hardly surprising 



3 




500 1000 1500 2000 2500 3000 



FIG. 1: 2D marginal 68% and 95% contours for z% and N v . 
Top: Linear prior for Zj. Bottom: Logarithmic prior, i.e., 
linear in log(l + Zi). 
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FIG. 2: 2D marginal 68% and 95% contours for z, and tlyih 2 , 
using the linear prior for z% and using N v as a fit parameter. 
Top: Full data set. Bottom: WMAP-5 data only. 



since N v ~ 3 allows the largest values of z% . 

In order to test more quantitatively how disfavored 
strongly coupled neutrinos are we have performed a high- 
precision run with N u = 3.046 and the more conservative 
linear Zi prior to calculate a sequence of progressively 
higher confidence limits. We find Zi < 1910 at 99% C.L. 
and 2230 at 99.7% C.L. At even higher confidence limits 
the Markov chains show signs of incomplete convergence 
and we refrain from quoting bounds. 

We also show 2D credible regions in the plane spanned 
by the matter density S1m = ^cdm + and in Fig. [2] 
where the conservative linear zi prior was used and N v 
kept as a fit parameter. In the top panel we have used 
the full data set as in Fig. [1] and find consistent results. 
In the bottom panel have used only WMAP-5 data and 
thus confirm with our method that WMAP-5 data alone 
do not significantly constrain Zj [28]. 

IV. LIMIT ON COUPLING STRENGTH 

Our approach of neutrinos recoupling at some red- 
shift Zi was motivated by a majoron-type interaction 
model where neutrinos interact with a new massless pseu- 
doscalar by virtue of a dimensionless Yukawa coupling g. 
In the framework of such a model we can translate our 
limit on z% into a limit on g in analogy to a previous pa- 
per by two of us [24j. When considering the scattering 
process the bound applies to any component of gij, the 
indices referring to the different neutrino flavors. The 



off-diagonal parts, however, are much more tightly con- 
strained by the decay process Vi — ► i/j<j) [24| . 

At z ~ 1500 the universe is matter dominated and to 
a good approximation H oc T 3 / 2 . Since for scattering 
the rate is Y ~ g 4 T, we can translate the condition for 
strong interaction, Y/H > 1, to a bound on g (24|. Since 
Y/H oc .g 4 ! 1 - 1 / 2 , and in the previous paper we effectively 
used Zi = 1088 to obtain g < 10~ 7 we now get g < 
lO-^lSOO/lOSS) 1 / 8 - 1.05 x 10~ 7 , i.e. a negligible 5% 
difference compared with our previous result. 

It should be noted that for masses below the recom- 
bination temperature, m < Tr ~ 0.3 eV, our bound 
applies equally well to the case of neutrino decay and in- 
verse decay [24j . In fact, we can now make the bound 
more quantitative by adding that at 95% C.L. the inter- 
action cannot be very strong before Zi ~ 1500. Roughly 
this translates to a bound on the lifetime (again scaling 
from our previous limit derived using z% — 1088) of 

r> 1.0 x 10 10 s (— ^— f. (2) 
V50meV/ v ; 

This limit is slightly weaker than before [24[ because of 
the slightly more conservative assumption about but 
it remains by far the most restrictive bound on invisible 
decays of low-mass neutrinos. 
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V. DISCUSSION 

We have updated bounds on the coupling strength be- 
tween neutrinos and a new, light pseudo-scalar, <p, using 
the latest cosmological data. Performing a slightly more 
refined calculation than in an earlier paper by two of 
us 



24j we find essentially unchanged constraints. 



One way to improve this limit in future is by actually 
detecting neutrino hot dark matter in cosmological pre- 
cision data. In decay scenarios involving massless pseu- 
doscalars and for a mass of 50 meV, the lifetime limit 
would improve by some six orders of magnitude [4^ |. 

Our more general conclusion is that neutrinos which 
are strongly interacting around recombination are 
strongly disfavored by data. The present data strongly 



support the conclusion that the cosmic neutrino back- 
ground (i) exists around the epoch of recombination and 
(ii) its fluctuations do have an anisotropic stress com- 
ponent. In future CMB data alone will likely suffice to 
reach the same or better sensitivity so that the bound on 
g can be expected to improve significantly [H, HH . 
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